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ABSTRACT 

We use the spectra of ~ 22, 000 early-type galaxies, selected from the Sloan Digital Sky Survey, to 
infer the ages, metallicities and star formation histories of these galaxies. We find clear evidence of 
"downsizing" , i.e. galaxies with larger velocity dispersion have older stellar populations. In particular, 
most early- type galaxies with velocity dispersion exceeding 200 km s" 1 formed more than 90% of their 
current stellar mass at redshift z > 2.5. Therefore, star formation was suppressed around this redshift. 
We also show that chemical enrichment was rapid, lasting 1-2 Gyr and find evidence that [Fe/H] is 
sub-solar. We study the robustness of these results by comparing three different approaches: using (i) 
Lick absorption line indices; (ii) fitting a single-burst stellar population model to the whole spectrum 
(lines+continuum); and (iii) reconstructing the star formation and metallicity histories in multiple age- 
bins, providing a method to measure mass-weighted ages and metallicities. We find good agreement 
between the luminosity- weighted ages and metallicities computed with these three methods. 

Subject headings: galaxies: general - galaxies: elliptical 



1. INTRODUCTION 

It is of general interest to theories of galaxy formation 
to determine observationally how early-type galaxies 
formed as this gives cl ues on how feedback was regulated 
in the s e systems (e.g. iBenson et all (I2002T) : iBaueh et~ 
1 2222 ) ; iHernaui st fc Sp ringe!! ll2003l l : iScannapieco fc Oh 
1 2004 i : ^Somerville et al J IgOOl : iNagamine et al. 
Il200l: iDe Lucia et all (|2006ft : iBower et al. I l|2006ft : 
ICroton et al. I |)200(tI) '1. The most direct route to do 
this is to observe the universe at different redshifts. If 
early- types are found up to a certain redshift, then one 
can determine how far back early-type galaxy formation 
occurred and what was the stellar mass and appearance 
of these galaxies. This indeed would be the best way to 
break the age-metallicity degeneracy, i.e. the property 
by which a galaxy can appear young or old due to 
complete lack of knowledge of its metallicity. However, 
as star formation becomes more intense in the past, a 
galaxy may hardly be classified as an early-type, thus 
making it difficult to find the progenitors of early-types 
by direct observation. Further, obtaining large samples 
of early-type galaxies at redshifts larger than one is very 
demanding observationally with the current generation 
of lOm-class telescopes. Nevertheless, there has been 
recent progress at finding the progenitors of nearby 
early-type gal a xies in the redshift range 1 < z < 2 (e.g. 
ICimatti et aH ll200l: IDaddi et all J200l: ILabbe et all 
(12005ft: IPapovich et all (120051): iCaputi et all 120061 ; 
iPapovich et alJ l)2006ft : iStern et all l|2006|) L however the 
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number of progenitors found is only about one hundred. 

One alternative approach to increase significantly 
the statistics, and be able to do environment studies 
ijSheth et al.ll200 6'l is to look at the local population of 
galaxies and try to infer from their stell ar populations 
the ag e s and metallic i ties o f the stars (e.g. iHeavens et alJ 
(j2004ft : |p"anter et alJ Q2006ft '). This approach suffers from 
two drawbacks: first, it can only tell us how the stars 
were formed, and not how the dark matter and stellar 
material were assembled. Second, recovering informa- 
tion from the integrated light of the stellar populations is 
subject to degeneracies: for certain features of the spec- 
tral energy distribution (SED) of galaxies it is not possi- 
ble to distinguish between an old stellar population with 
a certain m etallicity and a younger one with a higher 
metallicity ijWorthevI [19 94) . On the other hand, since 
galaxies are nearby they are not faint and therefore large 
numbers can be observed with small telescopes. Fur- 
ther, the use of the full spectral information significantly 
reduces, o r removes altogether, th e age-metallicity de- 
gencracy (IFanelli et al.lll987l IT992) IDunlop et aT]ll996T: 
Spinrad et alJl!997t IDornianetldTl2003l: IJimenezet alJ 
20041) . 

An interesting question is whether given enough spec- 
tral coverage and signal to noise in the SED, it is possi- 
ble to recover any evolutionary information from obser- 
vations of nearby early-type galaxies: the fossil record. 
This question has been ad dressed before based purely 
on theoretical grounds (e.g. IJimenez et alJ (|2004)). The 
question that we investigate in this paper is similar but 
we do it by using observations from the Sloan Digital Sky 
Survey (SDSS) and stellar population models that specif- 
ically include the treatment of non-solar scaled abun- 
dances. In particular, we want to know if the recon- 
structed star format ion history derived from MOPED 
ijHeavens et al.l 12000ft is consistent with the luminosity- 
weighted ages and metallicities obtained with other pro- 
cedures. Early-type galaxies are particularly suitable for 
this study since they have simpler star formation histo- 
ries than disk galaxies. 
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In a previous study ijHeavens et al.ll2004j) . we recov- 
ered the star-formation histories of a sample of SDSS 
galaxies, regardless of type/morphology. In this paper, 
we specifically analyze the subset of early-type galax- 
ies. Our main findings are as follows: the star forma- 
tion of early-type galaxies can be accurately recovered 
with MOPED, even if the number of parameters is sig- 
nificant. From the sample of early-type SDSS galaxies 
we find clear evidence for "downsizing", i.e. most mas- 
sive galaxies form their stars first. This effect is more 
pronounced for early-type galaxies than for the whole 
galaxy population. Furthermore, we find that more than 
60% of massive galaxies have formed 90% of their stars 
at z > 2.5. This places a strong constraint on when 
star formation has to be quenched in these systems. We 
show using hydro-dynamical models of early-type forma- 
tion how this can be achieved. Furthermore, we show 
how chemical enrichment took place in this systems. In 
total we use about 40, 000 early- type galaxies, a number 
that is many orders of magnitude larger than any direct 
observations at z ~ 2 — 3. 

The paper is organized as follows: §2 describes the 
sample used in this study. In §3 we describe the differ- 
ent methods used to determine the ages, star formation 
histories and metallicities of early-type galaxies. §4 and 
5 describe the comparison between the different meth- 
ods. We present the star formation history in §6 and the 
metallicity history in §7. We draw our conclusions in §8. 

2. THE SAMPLE 

We ha ve used the sample of early-type galaxies ana- 
lyzed bv iBernardi et all l)2006j) . The sample, extracted 
from the Sloan Digita l Sky Survey llYork et all |200C|) 
in its Data Release 2 ijAbazaiian et al.ll2005l) . contains 
over 40, 000 early-type galaxies selected for having ap- 
parent magnitude 14.5 < r < 17.75 with spectroscopic 
parameter eclass < 0, which gives a PCA component 
corresponding to no emission lines, and fracDev r > 0.8, 
which is a seeing-corrected indicator of morphology. In 
addition, only those objects with measured velocity dis- 
persions were selected, which translates in the spectra 
having S/N > 10 in the region 4200 - 5800 A. The sam- 
ple encompasses a redshift range 0.05 < z < 0.2, which 
corresponds to a maximum look-back time of 2 Gyr. 

3. THE METHOD 

3.1. Lick Indexes 

IBernardi et alJ 1120061) used Lick inde xes llWorthevI 
119941 iTraeer et al. 1 11991 iThomas et al~l l2003j) to de- 

rive the luminosity-weighted age and metallicity of the 
galaxies in the sample. In order to be able to ob- 
tain enough signal-to-noise in the spectrum to accu- 
rately measure the absorption features, they stacked 
the spectra that had similar properties. This reduced 
the sample to 925 composite spectra. IBernardi et al.l 
(2006) use the Mgb, <Fe>, H/3 and IPy F Balmer ab- 
sorption lines and a— enhan ced stellar population models 
l)Thomas e\ al. |f2QQ3j |2Q 04fl to derive ages and metallic- 
ities. See lBernardietaTTl|2006j) for more details. 

3.2. Full spectral fit 

Our second technique takes advantage of the large 
spectral range covered by the SDSS spectra, specifically 



the measured flux in the region 3500 - 4000 A. This re- 
gion is especially useful wh en trying to disentan gle the 
age- metallicity degeneracy ll Jimenez et alJl2003) . Fur- 
thermore, a fit to the full spectrum takes advantage of 
the sensitivity of the cont inuum to metal blanketing . We 
have used the models bv lBruzual fc Charlotl l|2003fl at a 
resolution of 3A. We then performed a \ 2 minimization 
to find the best fitting model in age and metallicity, so 
called single stellar populations (ssp). We note that we 
performed the fit to individual spectra, as we find that 
the average S/N of the spectra (above 10 per resolution 
element) is sufficient to yield a meaningful fit. 

3.3. MOPED fit 

The third and final meth od uses the MOPED algo- 
rithm l|Heavens et al.ll2000D that goes one step beyond 
ssp fitting and recovers in a minimally parametric way 
the star for mation and metallici ty histories of a given 
galaxy (e.g. iHeavens et all 120041) . The re cent DR3 re- 
lease f rom the SDSS has been analyzed bv lPanter et alJ 
(2006) using MOPED. It is from this study that we se- 
lect those galaxies that overlap with the Bernardi catalog 
of early-types. This results in about 22,000 galaxies in 
common. 

In brief, MOPED uses 11 bins, equally spaced logarith- 
mically in lookback time, to describe the star formation 
history of a galaxy. In each bin there are two unknowns: 
the strength of the burst and its metallicity, that are 
searched from a set of synthetic stellar population mod- 
els, to produce the best fit to the observed galaxy. In 
addition there is a further parameter that determines 
the amount of dust present in the galaxy toda y . For 
this study we have used the iBruzual fc Charlotl |2003) 
models at 3A resolution as the ssp describing the stellar 
population in each of the bins. 

One obvious concern with MOPED is that with a 
search in such large parameter space, degeneracies may 
be significant. Although this could be the case on an 
individual basis. IPanter et ail l)2003l) showed that this is 
not the case for stellar populations with old comp onents 
(age > 8 Gyr). Furthermore. IPanter et all (|2003L 12006]) 
demonstrated that no degeneracies remain for large sta- 
tistical samples. One clear advantage of MOPED is that 
it can provide the star formation history of the galaxy 
and therefore the possibility to calculate a mass- weighted 
age and metallicity as opposed to luminosity-weighted 
quantities as the other two methods. 

In the following sections we show that when one com- 
putes the luminosity-weighted age and metallicity from 
MOPED they agree well with those obtained from the 
Lick or full spectrum methods. This shows that de- 
spite the search in such a large parameter space (23 pa- 
rameters), there are no significant degeneracies in the 
MOPED solution. 

4. THE METALLICITIES OF EARLY- TYPE GALAXIES 

Fig. H shows the comparison between the luminosity- 
weighted meta ll icity d erived from the Lick indexes in 
IBernardi et"aH l)2006l) and from MOPED. The Lick 
m etallicities are com puted with the a— enhanced models 
of lMarastonl ((2005) while th e MOPED metallici t ies ar e 
for the solar-scaled models of Bruzu al &: Charlotl (J2003) . 
For the MOPED method we compute mass-weighted 
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Fig. 1. — Comparison between the luminosity- weighted metallic- 
ities obtained from the Lick indexes and the luminosity-weighted 
metallicity (see text) obtained from MOPED. There is good agree- 
ment between the two methods. 




Age SSP [Gyr] 

Fig. 2. — Comparison between the luminosity-weighted age ob- 
tained with fits to the full spectrum (Age SSP) and MOPED. There 
is a slight trend for MOPED ages being older by about 5% for the 
oldest ages. 



metallicities by adding up the 11 bins weighting by the 
luminosity of each component in the SDSS band r. 

In order to convert the solar-scaled met allicities to the 
a— enh anced scale, we use the formula bv lTraeer et al.1 
(120001) . To do the com parison with the composites from 
iBernardi et alJ l)2006|) we have averaged the metallicities 
of the MOPED spectra that correspond to a given com- 
posite in the Bernardi et al. sample. 

We find good agreement between the two methods. 
There is a slight systematic shift at high-metallicities 
of about 0.1 dex. The MOPED corrected metallicities 
tend to be lower than the Lick derived ones. The la dis- 
persion around the one-to-one line is of about 0.05 dex. 
The metallicities for all galaxies are super-solar, with no 
galaxies below the solar value. The maximum metallicity 
is about twice the solar value. 
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Fig. 3. — Fraction of galaxies in the SDSS sample as a function 
of velocity dispersion that have formed more than 90% of their 
stars at a look-back time greater than 11 Gyr and 8.5 Gyr (solid 
lines) . The dashed lines are the same but now for galaxies that have 
formed more than 70% of their stars. Note the strong correlation 
between fraction of galaxies and velocity dispersion which is clear 
evidence of "downsizing" for the early- type population. Also note 
that for log(av /kms - 1 ) > 2.3 more than 50% of the galaxies have 
formed 90% of their stars at z > 2.5. This fraction increases to 
90% of galaxies that have formed more than 90% of their stars 
at 2 > 1.2. This indicates that a significant fraction of massive 
early-types have shut-off stars formation at z ~ 2.5. 



5. THE AGES OF EARLY- TYPE GALAXIES 

Fig. |21 shows a comparison between the luminosity- 
weighted ages obtained from the single fit to the whole 
spectrum and the MOPED luminosity-weighted ages (i.e. 
summing-up the ages of all the 11 bins weighted by their 
luminosity in the SDSS band r). 

This is the more direct method to determine the extent 
degeneracies can play in the MOPED solution. Because 
the two methods use the same se t of stellar populations 
models ijBruzual fc Charlotl2003|) . the comparison is free 
of systematics. It gives a measurement of MOPED de- 
generacies because it compares a fit with just two pa- 
rameters, age and metallicity, with the 23 parameters of 
MOPED. 

The agreement is reasonable. There is a tendency for 
the MOPED ages greater than 10 Gyr to be older than 
the SSP ages by about one Gyr. However, this can be un- 
derstood from the fact that the last bin of age in MOPED 
is quite wide (10 to 13 Gyr) and therefore the uncertainty 
is of the order of 20%. As Fig.|2]shows, and also discussed 
in § 6 below, the weighted age is dominated by this last 
bin, so the intrinsic uncertainty in the MOPED bin dom- 
inates the error in the comparison. In the comparison we 
chose 11.5 Gyr as the 'age' of the bin. 

The comparison with the Lick indexes ages is not 
as straightforward, because the Lick index ages were 
calculated with a— enhanced models. Fully-tested 
a— enhanced models that can be used to fit the full spec- 
trum (lines+continuum) are not yet available. 

6. THE STAR FORMATION HISTORIES OF 

EARLY-TYPE GALAXIES 

We can now look in detail at the star formation his- 
tories provided by MOPED. In particular, we are inter- 
ested to know what the typical star formation history of 
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an early-type is and how it relates to the mass of the 
galaxy. 

Fig. |3 shows what fraction of early-type galaxies in 
the SDSS have formed more than a given fraction of 
their current stellar mass above a given look-back time. 
The solid lines show the fraction of galaxies that have 
formed more than 90% of their current stellar mass at 
look-back times larger than 11 and 8.5 Gyr. First note 
the strong correlation with velocity dispersion which 
is a clear sign of "downsizing" for early- type galaxies, 
i.e. star formation o ccurs in more massive galaxies at 
earlie r times (see e.g. ICowie et a l.lll999l iBernardi et al.l 
2003 IHeavens et all 12004 IKodama et al. I 120041 
Tho mas et al. 1 120051 iTanaka et al. 1 120051 llreu et al. 
20051 Ivan der Wei et al. I 120051 Duneau et al" 
20051 IBernardi et all 120061 iBundv et al. I [ 20061: 
di Serego Alighieri et al. I 120051: iCimatti et all 120061: 
Neistein et al.ll2006HMouri fe Taniguchil l2006'). Also, for 
galaxies with log(av/kms -1 ) > 2.3 more than 50% of 
the SDSS early-types have formed more than 90% of 
their stars at z > 2.5. This fraction increases to 80% at 
z > 1.2. Clearly large galaxies {ay > 200 km s _1 ) are 
shutting off their star formation at very early epochs. 
The dashed lines show a similar trend for a fraction 
in stars larger than 70%. For the whole population 
of massive early-types, there is very little stellar mass 
being added below z — 1.2. 

Therefore most of the star formation in massive early- 
types took place at z > 2 and in a sh ort burst of du- 
ration 1 — 2 Gyr llimenez et al.l ll999l h The hydrody- 
namical models in I Jimenez et all l)1999|) also explain the 
shut-off in star formation after just 1 — 2 Gyr as SN feed- 
back that heats the gas beyond the escape velocity of 
the dark matter halo. These models do not include any 
AGN feedback yet seem to reproduce well the trends for 
early- types in the SDSS. The clue to their success is the 
early assembly of the dark matte r halo . Similar trends 
were a l so found inlBlanton et al.l (|2000f) : iGranato et al.l 
(200i|): ljmTenez et all l|2005|) . 

Because the MOPED star formation history points to 
most of the star formation in early-types taking place at 
t > 11 Gyr, it is not surprising therefore that luminosity- 
weighted ages and m ass-weighted ages are very similar 
( Halman et al.ll2006ft . 

The MOPED ages are calculated over bin ages that 
are somewhat large (about 20% of the age) and they are 
not accurate enough to see the evolution of the red en- 
velope in detail as a function of redshift. However, it 
is interesting to look at the full-spectrum fit luminosity- 
weighted ages as a function of redshift. This assumes 
that the star formation took place in a single burst, but 
as we have seen above this is suggested by the MOPED 
stars formation history. Fig. 0] shows the age of the old- 
est early- type galaxies as a function of observed redshift. 
The plot shows a contour where the white region includes 
99% of the galaxies. For reference, the solid line is the 
age-re dshift relation for th e currently favored LCDM uni- 
verse llSpergel et al.ll2006l) . It is left for future work to 
provide a robust dete rmination from the se data of the 
Hubble constant (e.g. iSimon et all l)2005|) 1. 

7. THE METALLICITY HISTORIES OF EARLY-TYPE 
GALAXIES 
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Fig. 4. — The age-redshift distribution for early- types in the 
SDSS. The white region contains 99% of galaxies. There is a clear 
old edge: higher-redshift galaxies are younger than lower-redshift 
ones. The solid line corresponds to the age-redshift relation for the 
LCDM universe. 
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Fig. 5. — Mass- weighted Fe abundance as a function of velocity 
dispersion for different look-back times. Note the trend of increas- 
ing metallicity with increasing velocity dispersion at all epochs. 
Furthermore, there is an increase in metallicity from t > 11 Gyr 
to t > 8.5 Gyr and then the metallicity remains almost constant, 
consistent with a picture in which there is n o further inflow of fresh 
gas. Also, we show the Fe value derived by Bernardi et al. (2006) 
from their Lick-index study at t = Gyr. 



We now turn our attention to the evolution of metal- 
licity in early-types. Fig. [3] shows the evolution of the 
mass-weighted metallicity, computed from the MOPED 
fit, for the SDSS early-type sample as a function of ve- 
locity dispersion and look-back time. As for the mass- 
weighted age, there is a correlation between metallicity 
and velocity dispersion. The more massive galaxies are 
not only the ones that contain the oldest stars but also 
the more metal rich. 

This is not too surprisi ng. In fact, a simplistic model 
(e.g. Uimenez et alJ l)1999l) ) in which one assigns a deeper 
dark matter potential to the most massive galaxies will 
predict just that, the more massive objects are able to 
trap the ISM gas longer due to the higher escape veloc- 
ity from the dark halo. This in turn means that gas in 
the more massive halos is more enriched than in the less 
massive ones. 

Furthermore, this model predicts that star formation 
will stop after 1 — 2 Gyr due to the gas been heated-up 
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above the escape velocity by SN explosions. From Fig. |3] 
we can see that this is the case, metal enrichment finishes 
at a look-back time of 8.5 Gyr (z > 1.2). Another inter- 
esting fact is that the average [Fe/H] of the population is 
slightly below solar when mass weighted. Note that the 
overall metallicity is however super-solar (Fig. GJ. This 
is also what one expects if star formation is quenched in 
1 — 2 Gyr. In this case the Type-la SN are not yet formed 
however they are the peak producers of Fe. 

8. CONCLUSIONS 

We have presented further evidence of clear 'down- 
sizing' for the SDSS early-type galaxies. We have also 
shown that these galaxies have formed a significant num- 
ber of their stars at z > 2.5, especially the most mas- 
sive ones. Furthermore, we have presented evidence that 
the build-up of stellar metallicity in early-type galax- 
ies monotonically increases with time and star forma- 
tion has been quenched just after 1 — 2 Gyr. This can 
be explained by a simple model in which SN heating 
of the interstellar medium is sufficient to suppress fur- 
ther star formation by expelling it from th e dark mat- 
ter halo of the galaxy. Uimenez et alJ 1)20051) noted that 
in the LCDM cosmology these massive halos are abun- 
dant at z ~ 3. They constitute the most natural can- 
didates to host the progenitors of the SDSS early-types. 
The problem that remains to be solved in simulations 
is to suppress star formation by subsequent gas-rich 
merge rs. This is precisely what A GN feedback provide s 
(e.g. iScannapieco & Ohl (I2004T): iBqwer et al. I (|2006j) : 
iCroton et al. I (|2006|) ; lDe Lucia et all 1)2006^ . 

Also, we have studied in detail different methods to 
recover the age(s) and metallicity(s) for the integrated 
stellar population of the SDSS early-types. We have first 
compared the two most commonly-used methods to re- 
cover the luminosity-weighted age and metallicity: fitting 
the whole spectrum with a single-burst model, and us- 
ing only Lick absorption line indices. We have shown 
that the two methods provide similar answers. Further- 
more, we have then compared these two methods to the 
MOPED algorithm, which does recover a (mostly) non- 
parametric star formation and metallicity history of the 
galaxy. Despite the larger freedom in MOPED, we have 
shown that the MOPED luminosity-weighted age and 
metallicity agree well with the above methods. Thus de- 
spite the larger number of variables MOPED uses, there 
does not seem to be significant degeneracies in its solu- 
tion. 

The large number of early-types discovered by the 
SDSS and the fossil record analysis provide a valuable 



window into the formation of this objects at z ~ 2 — 3 
when the universe was only 1/5 of its current age. We 
have been able to use more than 40, 000 early-type galax- 
ies for our study. Such a large number of galaxies is 
not currently available at z ~ 3, maybe we know of a 
few handful of secure early-type progenitors by direct 
observation. As an example of the usefulnes s of such a 
large sample, in a companion paper ijHaiman et al.l2006j) 
we use this fossil sample to constrain the evolution with 
redshift of the physical parameters that determine the 
growth of massive black holes at the centers of galaxies. 
With the tight constraints imposed by the fossil record, 
it is clear that a significant number of EROs should be 
detectable at such redshifts with the next generation of 
telescopes. 
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